


Institutional Archive of the Naval Postgraduate School 





Calhoun: The NPS Institutional Archive 


DSpace Repository 


Theses and Dissertations 


l. Thesis and Dissertation Collection, all items 


1972 


Cool-down problem of a solid propellant 
rocket motor incorporating dewetting. 


Wellmann, Donald Albert. 


Monterey, California. Naval Postgraduate School 


http://hdl.handle.net/10945/16056 


Downloaded from NPS Archive: Calhoun 


MUS | DUDLEY 


WN] | orar 


http://www.nps.edu/library 





Calhoun is the Naval Postgraduate School's public access digital repository for 

research materials and institutional publications created by the NPS community. 

Calhoun is named for Professor of Mathematics Guy K. Calhoun, NPS's first 
appointed — and published — scholarly author. 


Dudley Knox Library / Naval Postgraduate School 
411 Dyer Road / 1 University Circle 
Monterey, California USA 93943 


COOL-DOWN PROBLEM OF A SOLID PROPELLANT 
ROCKET MOTOR INCORPORATING DEWETTING 


ТҮТІН Albert Wel mann 








NAVAL POSTGRADUATE SCHOOL 


Monterey, California 





THESIS 


Cool-down Problem of a Solid Propellant 
Rocket Motor Incorporating Dewetting 


by 


Donald Albert Wellmann 





Thesis Advisor: | G. H. Lindsev 


June 1972 


Anproved бол public release; distribution unlimited. 


1147655 





Cool-down Problem of a Solid Propellant 
Rocket Motor Incorporating Dewetting 


by 


Donald Albert Wellmann 
Lieutenant, United States Navy 
B.S., United States Naval Academy, 1964 
M.S., Naval Postgraduate School, 1971 


Submitted in partial fulfillment of the 
requirements for the degree of 


AERONAUTICAL ENGINEER 


from the 


NAVAL POSTGRADUATE SCHOOL 
June 1972 





ABSTRACT 


This thesis undertook a macroscopic, theoretical analysis of a long, 
circular, solid propellant, rocket motor undergoing thermal loading. It 
was shown that classical plasticity's von Mises yield condition adequately 
described dewetting of solid propellant which was modeled to have linear 
strain hardening characteristics. The mechanical properties of the 
propellant permitted a closed-form solution for the stresses and the 


strains resulting from cooling the motor from 160°F to -60°F. 
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I. INTRODUCTION 


The problem analyzed in this thesis was that of a long, cylindrical, 
solid propellant rocket motor (see Figure 1) containing a circular bore 
and with a thin steel case bonded to its outer radius. The motor was 
thermally loaded by cooling it from 160°F to -60°F, which represented 
the extremes that might realistically be expected during its lifetime, 
from cure at the upper temperature to storage at the lower one. The 
cool-down rate was assumed to be slow enough to maintain equilibrium 
conditions throughout the propellant; i. e., temperature did not vary with 
radius. The stresses in the motor resulted from the difference in the 
coefficients of thermal expansion between the case and the propellant. 
Because the propellant had the higher of the two coefficients, it tended 
to contract more than the case; thus an applied load was induced at the 
outer boundary of the propellant. 

The propellant was considered to be an isotropic, solid material; 
although in reality it was a composite of solid oxidizer particles 
imbedded in and bonded to a carbon fuel matrix. After an initial linear 
region, nonlinear stress-strain behavior under loading resulted from the 
progressive failure of these bonds at the interfaces between the fuel 
matrix and the oxidizer particles. Voids or vacuoles were thus created 
throughout the material, and their volume increased with further loading. 
This phenomenon is called dewetting. 

Dewetting causes nonlinear material behavior and has many similari- 


ties with yielding in classical plasticity theory, as was pointed out 








by Lindsey and co-workers, (4) who developed an isotropic theory for 
dewettable solids. In it, the strain increment was decomposed into an 
elastic component and a dewetted component. 
~ 4 
€ E = 
d£; E dE: sh d£ (1) 
They further developed from Drucker's postulate that the dewetted strain 
increment was normal to the dewetting surface. 
dei; zu DECS; ,E:) (2) 
O 6 
where f=o describes the dewetting surface. From this theory, Wood (7) 


showed that the general form of f was 


и for I 20 


J,- B 
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where Jy was the second deviatoric stress invariant, I, was the first 


stress invariant, and A and B were mechanical properties. Wood experi- 
mentally measured these parameters for a PBAN propellant and found that 
A = -0.0026 and B = 180. 

Lindsey and Wood (4) described the dewetting surface as a circular 
cylinder about the hydrostatic axis in principle stress space. The 
positive end of this cylinder was capped by a body of favole Ion (see 
Figure 2). Using a space defined by an axis system consisting of L> 
which was coincident with the hydrostatic axis and J, which was perpen- 
dicular to the hydrostatic axis, Wood's dewetting criterion became 
intersecting straight lines with an extremely small negative slope of 
-0.0026 and an intersection with the J, axis of 180 , as shown in 
Figure 3. Taking the elastic propellant properties of Wood's material and 
the. elastic solution for an encased cylinder subjected to thėrmal cool- 


down, the intersection between the thermally induced stress field and 


the dewetting surface was determined by calculating the ratio of Jo and 








and г, and finding its intersection with the dewetting surface, 


(See Figure 1) It was found that only 1 1/3% error resulted in assuming 
the dewetted surface was a cylinder for the thermal problem. This was 
quite acceptable, especially in view of the great mathematical simpli- 
fication which resulted. Thus, Wood's dewetting criterion for this 
problem was approximated to be 

J + B(T) (4) 
where for the thermal problem, B is a function of temperature. It was 
observed that this dewetting criterion was identical to the von Mises 
yield criterion of classical plasticity. 

The goal of this thesis, then was to perform a stress and strain 
analysis using Wood's dewetting criterion as simplified accounting for 
the nonlinear behavior caused by dewetting, on a long, circular, solid 
propellant rocket motor. 

Prior theory had considered propellant to be completely elastic; 
however, the phenomenon of dewetting was expected to occur at ‘the 
inner bore initially and to expand outward as loading was increased. 
After yielding occurs in either dewetting or plasticity theory, the 
linear stress-strain law is no longer valid and must be replaced by а 
flow rule. Both theories insist upon a similar form of the flow rule 
and for the dewetting surface of equation (4), the flow law of equation 
(2) becomes 

dE, = Su (5) 
where dE; is the increment of dewetted strain, Se; is the deviatoric 
stress component, andA is a scalar EES the loading history. 
Mendelson (5) shows how À can be written in terms of Se an equivalent 


stress, and d£, an equivalent dewetting strain increment: 


A = 3 d& (6) 
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In this analysis both the dewetting function and the propellant's 
mechanical properties (except for Poisson ratio which was assumed 
constant) were functions of temperature. These properties were determined 
experimentally for a PBAN propellant manufactured by United Technology 
Center, Sunnyvale, California. Thus in the solution as the temperature 
loading was changed, it was possible to account for variable dewetting 


requirements as well as variable stress-strain relationships. 
% 








II. BACKGROUND 


Because of the close analogy between dewetting and plasticity, the 
pertinent literature of plasticity was searched for solutions to cylinder 
problems. It was found that the subject of partially plastic, thick- 
walled tubes received much attention from many authors during the late 
forties and the fifties, In general, the vast majority of the work 
dealt with internally pressurized cylinders, which were then of current 
interest and were more amenable to solution than those problems which 
deal with thermal loading. 

In 1957, Steele (6).made a review of several theories for partially 
plastic thick-walled tubes. He stated: 

"The mathematician's approach to the problem under 
consideration is fundamentally different from that of 

the engineer. The former is interested in an accurate 

knowledge of the stresses and strains in terms of the 

applied load or loads for all points in an idealized 

body. In general, the engineer seeks an estimation of 

the variation of factors with applied loads, which 

are of significance only in the functional design of 

the part.... In the case of the thick-walled cylinder, 

interest lies in the deflections of the bore and the 

outer surface for specific internal pressures." 

As an example, Hill, Lee and Tupper (2), in a rigorous work, showed 
that Nadai's approximate solution for the axial stress was in error by 
sixty percent, although the deflectionsat the bore and outer surface 
correlated closely. Rigorous solutions such as Hill's, Lee's and 
Tupper's, did not yield a closed form solution; and, although mathema- 


tically correct, they were so complicated that their usefulness to the 


engineer was minimal. 








A general outline of the relationships required by the thick-walled 
cylinder problem is given below. These are required for any mode of 
loading, and apply to elastic, plastic, or dewetting theories equally. 

(1) Radial Equilibrium 

(2) Strain Compatability 

(3) Stress-strain Relation 


(4) Continuity of radial stresses and displacements across the yield 
surface 


(5) Boundary Conditions 

(6) Compressibility 

@) yield Criteria 

(8) Yield zone stress-strain relations 

The boundary conditions are fixed by the geometry and the physical 
environment. It was assumed that at the inner bore the radial stress 
was zero and that the outer surface had an induced pressure caused by 
the restraining effect that the case applied to the propellant as the 
motor was cooled. Plane strain was the assumed end condition. 

For mathematical simplicity in many plasticity ns, the material 
has been assumed to have the same compressibility characteristics in the 
yielded region as in the elastic region. Furthermore, the material 
has been taken to be incompressible, which for metals can introduce 
considerable error. However, this assumption is quite valid when dealing 
with solid propellant, as Wood's experimental evidence has shown. Thus, 
it was assumed that the material was incompressible in both the dewetted 
and the undewetted regions. This gave considerable mathematical 
simplification in that Poisson's ratio, JJ, was assumed to be a constant 
and equal to one half and that the sum of the dewetted strain increments 


was equal to zero. 





i 3 i 

dta dta + аё, ZU (7) 
This assumption permitted the simplification of the yield criterion 
and the stress-strain laws to useable forms. i 

Two yield criteria were frequently encountered. The von Mises, 

where 
K 
J, ^ eum (8) 
and the Tresca yield condition, in which 
|6- 94 = (9) 

The Tresca condition offered great mathematical simplicity for the 
cylinder problem facilitating solutions that were otherwise out of reach, 
provided that the relative magnitudes of the three principle stresses 
were known a priori. It was shown that in this problem, the axial stress 
was the middle stress, thus there was no need to calculate it. Addi- 
tionally, its contribution to plastic yielding was zero. 

On the other hand, the von Mises criterion required no prior 
knowledge of the stresses, but it yielded nonlinear differential 
equations. Due to the fact that WE have demonstrated 
experimentally that the von Mises criterion provides a closer correla- 
tion with a description of yield surface than the Tresca condition, 
solutions have been sought and obtained for the re problem using 
the von Mises criterion. The thermal problem, however, was more difficult 
in that a thermal term was added to the stress-strain equation. The 
addition prevented a simple closed form solution using the von Mises 
criterion; however, it had no effect on the Tresca condition. Consequently, 


the Tresca condition, even though less accurate, was frequently chosen 
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over the von Mises criterion. Recall, though, that from dewetting 
theory Wood found the dewetting criterion to be of the von Mises form. 
It was therefore desired to obtain a solution using the von Mises 
surface. 

There are two plastic flow rules that have been used exclusively the 
Prandtl Reuss and the Hencky: 


Prandtl-Reuss 


i cox CD] 
e 
LET EPT s e ey (10) 
Ж 
NS 6% S (Set SN 


where dE» is the equivalent plastic strain and Se the equivalent stress. 


Hencky 
E x = § 6-4 Sor 6} 
о 2-і (4-24-64) (11) 
Е = 


z PS A + (e+ Sey} 

The Prandtl-Reuss equations deal with incremental theory where an 
increment of plastic strain is so small that it can be considered linear. 
The total plastic strain is achieved by integration. An expansion of 
equation (5) for dewetting materials shows it to be identical to the 
Prandtl-Reuss equations. 

Hencky postulated that the total plastic strain would have the same 
type of relationship as the incremental such that the present state of 
plastic strain was independent of its loading path. This is not true 
in general, and the Hencky relationships have fallen into disuse. 
However, under special circumstances of linear strain hardening and a 


monotomically increasing load, the two flow rules coincide (5). 


LI 





Steele presented a summary of the combinations of these parameters that 
have been tried by different researchers. It is repeated here in Table I. 
Note that most of the solutions required numerical analysis. Any 
numerical solution necessitates a loss of generality and makes a para- 
meter study more costly and difficult. Thus a closed form solution was 
desired. By making simplifying approximations which were the result 
of either the special geometry assumed or the unique mechanical properties 
of solid propellant, it was possible to obtain a closed. form solution 


to the thermal problem using a von Mises dewetting surface. 


12 





ІІІ. THEORETICAL ANALYSIS 


Within the above mentioned framework and background, the theoretical 
analysis was undertaken. Figure 2 describes the geometry of the problem, 


and the governing equations are given below. 


Equilibrium 


46: ль (12) 
ШЕ e | 
B r 
Strain Compatability 
dto өсі =0 (13) 
dr r 
Stress-Strain 
4 к a 
Ее = = 6 P(S2 EY + AT + Er (14) 


Eo = $ Í Sp P (Set EJ XATT Ea 
Ez > z f Sar let бәр өле E, 


е 
where о are the dewetted strain components. 
d 2 


The problem was well structured ín 
equations to accommodate the unknowns; 
derived. 


By applying the plane strain 


bility condition of the dewetted zone, 


that there were sufficient 
thus, a unique solution could be 
assumption and the incompressi- 


the last of equations (14) was 


solved for the axial stress 


Сар 


: 1/2(6,+6 ө)-Еот + E(& +E.) 


(15) 

Using this definition of the axial stress, the first two stress-strain 
relations (equations (14) were substituted into the strain compatability 
equation (equation (13)) yielding a differential equation in terms of 


So and бы. The tangential stress, Gp, was eliminated through the use 


of the equilibrium equation (equation (12)), yielding 


13 





(16) 


= fa ans теке). ; Se] à 3 (En еу) 

г : : 

a6: (£e -£). 0 
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Lo [no 


+19 
“та 


Integrating thís equation twice yielded an expression for the radial 


stress. 


Va 4 5 
6 = -EQ,, 2E SEE dr — Ef aredo HS (17) 
r v? 3 “ar mE 


where C1 and C) were constants of integration. Analytic integration is 
only possible when Poisson's ratio is one-half, which is the actual 
value for propellant. 

The temperature distribution was assumed to be constant as previously 
mentioned. Consequently, the temperature integral of equation (17) was 


readily evaluated. 


ee 
6 2. ЕС + AE E Es An _ Eex(- E r С, (18) 
r et 3 26 r 


The two constants of integration were evaluated by applying the boundary 
conditions: 
at r = a $ 


at r=b б. = р’ 


0 


г 


(19) 


where p' was the pressure at the boundary between the 
propellant and the case and was caused by the difference between the 
case's coefficient of thermal expansion and that. of the propellant. 
See Appendix A. Application of the first boundary condition yielded 
ge 
С, rO (20) 
Applying the second boundary condition at r = b, gives 


E RP, Ear P: 
c= EY |--p - q) tdo + ( bat )] (21) 


Substituting (19) and (21) back into the expression for 6, yielded 


¿EA [y 2 Eto dr E AS SD 


~ (par 
at ers - £ax(1- € 
= La 6) (22) 
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The need was now to solve the remaining dewetted strain integrals. 
As mentioned above, it had been assumed that the Lindsey dewetted 
stress strain equations applied and that solid propellant possessed 
linear strain hardening characteristics. Expanding equations (5) gives 
d£, - e^ $ = KERTA, 
diese cto f so- т (6,69) 
e 


As previously stated, for monotonically increasing loads and a linearly 


(23) 


strain hardening material, the Prandtl-Reuss equations degenerated to 
the Hencky equations. Within this approximation, therefore, it was 


assumed that the Lindsey equations took the form of 


E. Е a CERTA 


: t i 
Eg * p» KA Gea 
From those equations it is seen that 
^ gz3 f (6. 
En ea 2 Ea с c 2) (25) 
where e Y LEN rt (Er ^ 
Bg = (Ets СС. e Y ( с te 
|; (26) 
8 LL (BS (6-60) 
6. =43Ј, Е ¡ISS | 
The axial stress, ‚ was replaced by equation (15), and after expanding 


the terms and regrouping 
Ge [y [Ce sY- мала) | 
(27) 
Ап examination for order of magnitude of each term of the right side 
of equation (27) demonstrated that the square of the stress difference 
clearly dominated the entire expression. .The radial and tangential 


stresses were of opposite sign and such size that their difference was 


relatively large. The low modulus of the propellant relative to that 
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of some other material such as steel (see Figure 13) made the EAT term 
small compared to the stress difference. Additionally, from the elastic 
solution, it is seen that the radial strain is negative while the 
tangential En is positive. Moreover, the absolute magnitude of the 
radial strain is larger than that of the tangential strain; consequently 
the sum of the tangential and radial dewetted strain components must be 
negative, assuming that the direction of the dewetted strain increments 
are in the same direction as their elastic components. So it seems 
that the already smallÉ£AV term is reduced even more by the addition of 
this negative quantity. Thus, the second term inside the radial was 
neglected. The expression for the equivalent stress became 

Ge 7136, Sol (28) 
which simplification made possible a closed form solution without 
appreciable loss in accuracy. The terms neglected were estimated to 
constitute at most a 2% error. 

Using equation (28), equation (25) reduced to 

Ep- - = 2_% san (Sr So) (29) 
where Sgn Ke Se <O for the thermal problem. 

The significance of assuming linear strain hardening actually 
meant that the stress-strain curve was eaten by two straight 
lines, one having the slope of E and the second having the slope of 
mE as shown in Figure 4. From the figure an expression for the dewetted 
region was found to be 

= (1-m) Se- Ŝo) 
Ep = Com 5.25, | (30) 
By taking equation (28) and applying the equilibrium equation (12) it 


is seen that 
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Using equations (31, 30, 29) permitted direct integration of the 
dewetted strain integral in equation (22) such that the radial stress 


expression became 


E Bleed. én 
По en ae ЛЕККЕ 


+ 13 GEEN Gene | 


Two new unknown constants appeared in thís equation, Śr, and € 


The critical radius, г.» Was the radius of ene boundary between the 
dewetted and the elastic zones, and Ġ, was the radial stress at this 
boundary. Both of these were constant and were functions of the thermal 
loading only. G, was related to r, by setting r equal to r 


с іп equation 


(32). After substituting for p' 


LANE (BETA AJ Blu) l al — 
(Q-a) -Ab Бо C 4: 2 a 
ME ea ND m facem Gole lq 4 
c Um-3 Сс оу ŻOEN 
Te бас E e vn 


where p' and D came from the elastic solution. (See Appendix A.) 

It was noted that at the boundary between the dewetted and undewetted 
zones, the radial stress must be continuous across the boundary. 
Further, note at that location that the dewetted condition must also be 
met by the stresses in both the elastic and the dewetted regions. 

The yield condition, equation (1) became 

(6-- So) = 24 Bm- Say (34) 

The second term under the radical was considered negligible for reasons 
already given. The expression for Se and So were taken from Appendix A. 
The combined loading case was used, and the radius was set equal to the 


critical radius (тег), such that 


£T (d. à 


(Ge 
È (E-o 


(35) 
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Combining equations (34) and (35) and не 9 for ES yielded 


a BETA | e TEE -i) la 
шш NU ЛУ” (36) 
2 — (-u)bE | 
5 ЗЕКЕ 
It would have been possible to solve equation (34) and (35) for S. in 


terms of an unknovn r, and to set this new expression equal to equation 


c 
(33), solving for ro directly. However, such a course of action would 
have resulted in a transcendental, tenth order expression and would have 
been of little practical value. A simpler nest was to iterate with 
equations (36 and (33). Assume a value for 6, (say 6, = 0) and solve for 
ra bv using equation (36). This value for г. was used in equation (33) 


to solve for a better value of <, . This process was repeated until 
с 


the resulting values of ге converged. 


All unknown constants of:equation (32) were known and the complete 


expression for the radial stress was obtained: 


< Y^ (rab "Ferny ex i а 








КУ 


As 2 ! 

PL ос акан Беа) 
r (ba) Y e (s (38) 

Through the application of the equilibrium equation, the tangential 


stress was obtained directly. 


EI es decente] ned) 


(>% р Cie: aba” у р- Хз С (566: & p 


с (је о 


(39) 


The strains were obtained by application of the stress-strain 


equations. 
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IV. EXPERIMENTAL RESULTS 


There was a dual goal in the experimental work done. First, it 
was necessary to experimentally determine the constant, B(T) in equation 
(2), for the particular PBAN propellant used. This was accomplished by 
conducting uniaxial tensile tests over a range of temperatures in the 
gas dilatometer of United Technology Center, Sunnyvale, California. 
For uniaxial tension, Equation (2) reduced to: 

2 
Dem = So (41) 

The dilatometer (see Figures 5, 6, 7, 8) utilized the Perfect Cas 
Law as its basic operational principle. The specimen was mounted within 
a sealed chamber and pull-rods extended from either end of the sample 
through the ends of the chamber and attached to the cross-head of an 
Instron testing machine. Any change in the volume of the sample while 
under load resulted in a change in the chamber's pressure creating a 
pressure difference across a diaphram separating the test chamber from 
an adjoining reference cavity. A baratron pressure transducer converted 
this pressure difference into an electrical signal, which was fed to 
one pen of a two-pen X-Y plotter. The signal for the second pen came 
from a load cell in the cross-head of the Instron machine. The raw data 
thus obtained were in the form of two superimposed plots; one was the time 
history of the load, and the second was the time history of the dilata- 
tion. By using the constant cross-head speed and premeasuring both 


the cross sectional area and gage length, the data was reduced to plots 


of stress vs strain and dilatation vs strain. (See Figure 9) 
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Changes in the ambient temperature caused a problem since it 
affected the readings of the pressure transducer through the same Perfect 
Gas Law. In order to isolate the experiment from ambient temperature 
fluctuations and also provide a method of controlling the temperature 
at which the experiment was conducted, the dilatometer was enclosed in 
a temperature box, which possessed the capability of maintaining a 
preselected temperature to within a degree over the duration of the 
experiment. Although the ambient temperature within the temperature 
box was fairly constant, at subzero temperatures, heat leaked into the 
dilatometer probably by conduction through the push rods. A large 
temperature drift would develop after about one half hour and this 
dictated the maximum duration of any individual run. 

The second goal was to experimentally determine the mechanical 
properties of the subject propellant as functions of temperature; 
namely, the uniaxial dewetting stress, Ó, , and the elastic modules, E 
and m. (See Figures 10, 11, 12) 

The test samples were square with sides machined to approximately 
one half inch, and they had a gage length of approximately 2 3/4 inches. 
Redwood tabs were bonded to each end of the sample to make hardware 
attachment possible. Tests were conducted at 21?C, 3?C, -9?C, -21?C, 
and -52°C, with four tests being conducted at the three higher tempera- 
tures, two tests conducted at -21°C and one test conducted at -52?C. 
Experimental difficulties with the equipment at the subzero temperatures 
resulted in the reduced number of runs at the lower temperature. 


Experimental scatter of seven to ten percent was found in the dilatation 


and stress-strain curves. However, only about five percent scatter was 
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found among the parameters measured (е. , ЕЯ UM and &, ). Due to the 
fragile nature of solid propellant specimens, this was considered quite 
reasonable results. 

After reducing the raw data to the stress vs strain and dilatation 
vs strain plots mentioned above, the problem of identifying 6 developed. 
Three methods were available from the literature: (1) interpreting the 
stress which corresponds to the lift-off point of the dilatation curve 
as <, . This method provided results which were quite repeatable; 
however, the stress strain curve did not depart from its initial linear 
relationship at this point. (2) The dilatation curves eventually 
assume a constant slope. Physically the total amount of solid material 
which could dewet had — at this point and further dilatation 
resulted only from increasing the size of previously formed voids. The 
stress corresponding to the extrapolation of this constant slope line 
back to the zero stress point (pt. B), of Figure 13 would be considred as 
the yield stress, 6. . This could be physically interpreted as assuming 
that no dewetting occurred until the yield stress was reached; then all 
the dewetting occurred simultaneously. This was tantamount to assuming 
that the propellant was perfectly plastic. However, at the lower 
temperatures, specimen failure resulted before such a linear portion to 
the dilatation curve was reached. (3) The third method was to take the 
stress at the point where the stress-strain curve initially departed 
from a linear relationship. This method reflected the linear-strain 
hardening assumption, and provided consistent reproduceable results 
from sample to sample. It neglected the dilatation curve altogether. 


This last method was adopted primarily because it met the condition 
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of the linear strain hardening assumption. A comparison of these 
three methods was given in Table II, where it is seen that methods 
(2) and (3) yield similar results. 

The test speeds were chosen as slow as possible, yet fracture had to 
be reached within the time limit mentioned above. At the lower tempera- 
tures, it was assumed that the speeds used maintained the specimens in 
equilibrium and that the load did not produce viscoelastic effects. 

In summary, the yield stress, on. was taken as the stress at which 
the average stress Strain curve first departed from linearity for all 
runs at any one temperature. The modulus, E, was the average of the 
elastic slopes of the stress-strain diagrams for any given temperature. 
The strain hardening exponent, m, was the ratio of the slope for the 
least squares straight line fit of the non-linear portion of the stress 


strain curve, to the initial modulus. 
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V. RESULTS AND CONCLUSIONS 


A: RESULTS 

In order to compare the results predicted by dewetting theory with 
those obtained from linear elastic theory, a sample problem was worked. 
A PBAN propellant of current interest was used, and its mechanical 
properties were experimentally determined, (Figures 11, 12). The other 
required parameters were taken from Appendix II of Reference (4). 


For the propellant: P= 1/2 A= 2.21х1072 іп. /іһ./°Е 

pomeete о себштсазе венка .65x107?  in./inc./?F 
E, - 3x10^ psi 

Geometrical parameters: a= 3in. b= 12in. h=0.015in. 

The rocket motor was cooled from 160°F to -60°F, and both elastic 
and dewetted stress-strain analyses were made. Linear strain hardening 
was assumed which enabled oa ee to be performed in two steps: 
one elastic step and one dewetted step. The temperature was not 


incrementally stepped down. Equations (32) and (36) were iterated in 


order to find the critical radius, ro and the interface stress, Gr. 


These results were plotted in Figure (14). It was interesting to note 
that the growth of the yield surface was linear with temperature. 
Furthermore, dewetting did not begin until the temperature had been 
decreased to 10°F, and at the lower temperature limit the yield surface 
had penetrated to a radius of only 9 1/2 inches. The interface pressure 
also increased linearly as the temperature was lowered reflecting the 
influence of temperature on both the yield criterion and the propellant's 


mechanical properties. 
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The greatest effect of dewetting was upon the tangential and radial 
stresses. (Figures 15, 16) The tangential stress at the bore, for both 
the elastic and the dewetting theories, was plotted against temperature. 
It was seen that dewetting substantially relieved the stress. The 
radial stress was plotted across the motor's thickness for a characteris- 
tic temperature of -8.0°F. Again it was shown that stress was relieved 
within the dewetted region. The radial stress then approached the 
elastic solution as one moved deeper into the and eenen elastic 
material and away from the critical radius. However, a pronounced 
effect of the induced pressure at the inner radius of the undewetted 
region was still evident. 

Initially, the behavior of the tangential strain at the bore was 
quite surprising. (See Figure (17)) It had been anticipated that this 
strain would have been negative and would have increased in magnitude 
as the temperature was lowered. However upon looking at the elastic 
tangential strain equation (A-1), one sees that it consists of two terms. 
The first term represents the mechanical strain imposed by the case 
upon the propellant. It is positive and dominates the second term, 
the thermal strain, which is negative in the cool-down problem. However, 
as the temperature is lowered the mechanical properties of the pro- 
pellant change significantly, and the negative T term grows in relative 
magnitude, actually forcing the total strain to decrease. Thus, it was 
observed that the maximum strain occurred at about 40°F instead of at 


the lower temperature extremum. 
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B. CONCLUSIONS 

A point worthy of note is the curve fit of the m curve coefficient 
of linear strain hardening in figure 12. Equipment limitations at the 
subzero temperatures, which restricted multiple runs decreased confidence 
in the values of m. Consequently, the specific mechanical properties 
measured at -5.8?F and -60?F were used to calculate the bore hoop 
strain at those two temperatures, and they appear on the plot in Figure 14 
which uses the curve fit for m. It was concluded that a curve fit which 
averaged these E icc provided, a smoother extension of the elastic 
solution than a curve which connected each data point. Further tests 
would prove useful to more definitely fix the m curve. 

From a comparison of the elastic and the dewetting solutions for 
both stresses and strains it was concluded that dewetting has little 
effect on the strains. In an incompressible material the strains are 
fixed by geometrical constraints and not by mechanical properties. 
However, stresses are tied to the mechanical properties and thus are 
greatly affected by dewetting. 

It was EE concluded that a thermo-dewetting solution could 
be obtained in a relatively simple form and that better results for stress 


would be obtained by using dewetting theory vice elastic theory. 
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TABLE II 


TABULATION OF EXPERIMENTAL DATA 


Temp ЗА Е m 
2d" 94 806 220177 
3.5 118 983 2505 
-9 137 1600 „59 
zn 159 8840 „369 
5 231 23100 „539 
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TABLE III. 


COMPARISON OF METHODS OF READING 
TENSILE YIELD STRESS 


TEMPERATURE 5 (psi) 

(96) І II III 
Dod 46 84.5 94 
345 72795 107 118 

-9 105 ж 13% 

-2u 137 ж 159 

25265 201 ж 2З 


*Dilatation Curve did not reach constant slope. 


I. % was stress corresponding to the lift off point 
of the Dilatation Curve. 


IL. 6, was stress corresponding to extrapolation of 
of the constant Dilatation Curve back to the strain 
axis. 

III. ó. was stress at which stress-strain curve initially 


departed from a linear relationship. 
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Circular 
Cross-section 


Оз 


Dewetted Yield Surface in Principle Stress Space 


Figure 2, 
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Intersection of Thermal Stress Load 
with Dewetting Surface 


Figure 3, 
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Linear Strain Hardening 


Figure +. 
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APPENDIX A 
ELASTIC SOLUTION 


Williams, Blatz and Schapery (7) of the Guggenheim Aeronautical 
Laboratory (California Institute of Technology) analyzed solid pro- 
pellant rocket motors of the same geometry under many different loading 
modes. Their basic assumptions were that | 

l. the И pellant had material properties which were 

a. Isotropic 
b. Homogeneous 
d. Continuous 

2. the strains were small. 

Their solution to the thermally loaded, cylindrical, rocket motor, 


assuming a conditionwof plane strain, was: 





2 b 
E. WAS e), де(-%) (rede E e jede 
Е E оху U TIBDA, 
РЕ y 
Sr = - bp (V i ), eie S) MN + Inde E — (A-1) 
(№? -ох\ (lli -a\ м 2 ©. S n 





л о? b 
- wrk. (De EA), rs $ (whe (red, 


Z E 50 
AR DIS ТТ 
where 2 (1455 b 
EA —— т ATI IE 
4 = Ada о 
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D= On (Landa Live Bis Б 
b’-ar - h Ee 
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And for the case of an internally pressurized rocket motor their 
solution was: 
Bet t 
E. о? р: рг), рге - р. 
S (0 о b'-al- 
217 : ы? 2 те 
| a eer y кесе. (А-2) 
Са b'-al 
BC 


Pe = zer БИО i 
| D 


where = internal pressure 


In considering the effects of dewetting, the inner boundary condition 
changed. The inner radius of the elastic region was the same as the 
critical radius, i.e., it was equal to the radius of the dewetting 
surface. As the dewetting surface moved outward, the inner radius 
of the elastic zone also increased; furthermore the inner boundary was 
no longer stress free. From continuity of the radial stress across the 
dewetting surface, it was seen that the inner boundary now had a stress 
equal toÓr, induced upon it. This elastic solution was obtained by com- 
bining A-1 and A-2, and setting p, equal to Sr., and the inner radius, 
a, equal to the critical radius, r, 

( 2 
Ġr € (r- b tal = ES = ЗВ (< el ET(d- AJ 
А Е с Ууу 2D 


j O 


3g (rdv ЗЫ (пч) 
¿AA CO ETL Aah 
Gg Xe (86-2 (Y - Y VD 280 2 (A-3) 


Sg 


P? (i-e 
3 3 nó. 
РЕ z EAT ere D-Ż v. (UN LE 
c D 2. (2 v?y ht. 
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